Transport and scattering in inhomogeneous quantum wires 
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We consider scattering and transport in interacting quantum wires that are connected to leads. 
Such a setup can be represented by a minimal model of interacting fermions with sudden changes in 
interaction strength and/or velocity. The inhomogeneities generally cause relevant backscattering, 
so it is a prion unclear if perfect ballistic transport is possible in the low temperature limit. We 
demonstrate that a conducting fixed point surprisingly exists even for large abrupt changes, which in 
the considered model corresponds to a velocity matching condition. The general position dependent 
Green's function is calculated in the presence of a sudden change, and is confirmed numerically 
with high accuracy. The exact form of the interference pattern in the form of density oscillations 
around inhomogeneities can be used to estimate the effective strength of local backscattering sources, 
offering a route to design experiments where the effects of the contacts are minimized. 

PACS numbers: 73.63.Nm, 71.10.Pm, 73.40.-c 
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The description of transport through quantum wires 
has become a very well studied research area, as it di- 
rectly ties together conductivity experiments* ~— with cen- 
tral theoretical models in one dimension, such as the 
Landaucr formalism^ and Luttinger liquids £ Landaucr 
showed that the conductance of a single spinless non- 
interacting quantum channel is always finite and given 
by G = (1 — R 2 )e 2 /h, where R is the backscattering 
amplitude^ However, interaction effects change this pic- 
ture dramatically as R becomes temperature dependent 
through renormalization. In the pioneering work of Kane 
and Fisher it was found that a local perturbation in a 
repulsive Luttinger liquid is relevant, which results in 
a characteristic powerlaw dependence R oc T 9 " 1 where 
g < 1 is the interaction parameter^ A number of impu- 
rity models in one dimension have since been analyzed in 
detail 8 - - — and confirm that generically a local perturba- 
tion cuts the transport at low temperatures, unless the 
relevant operator is forbidden by symmetry. 

Nearly perfect connections between leads and wires can 
be achieved for different types of quantum wires^— which 
show quantized conductance at moderate temperatures. 
At lower temperatures, however, scattering in the connec- 
tions to the leads plays an increasingly important role. 
The higher dimensional contacts are effectively weakly 
interacting, so there has been great interest in describing 
the transport through a quantum wire attached to non- 
interacting leadsJ£~— The minimal model for this setup 
is a single channel of interacting spinless fermions, where 
the interaction parameter changes along the wire. Under 
the assumption that it is possible to use a hydrodynamic 
Luttinger liquid description for this model^ it would be 
expected that backscattering follows a non-trivial renor- 
malization behavior which is position dependent^ Even 
if the connections are effectively free from imperfections 
with a homogeneous lattice structure,— there have to be 
small regions of the wire, the junctions, where the in- 
teraction changes, which does not necessarily occur adi- 



abatically and will induce intrinsic backscattering. This 
immediately invites the question if it is ever possible to 
create a perfect connection in the low temperature limit. 
Indeed it is in general unclear how large such intrinsic 
backscattering is for generic junctions and if it is even 
justified to use an inhomogeneous Luttinger liquid de- 
scription in the first place for this setup, since the usual 
bosonization procedure assumes a translational and scale 
invariant theory. 

In this paper we consider the intrinsic backscattering 
in an inhomogeneous wire and show that a perfectly con- 
ducting fixed point can indeed be found by adjusting pa- 
rameters such as the velocities on both sides. At half fill- 
ing it is also possible to estimate the bare strength of the 
backscattering from the change in the local velocity field. 
The full correlation function of an inhomogeneous Lut- 
tinger liquid is calculated, which agrees with numerical 
simulations to high accuracy. This confirms the validity 
of the inhomogeneous theory and proves that a low en- 
ergy conducting fixed point can exist even in the presence 
of large abrupt jumps. 

A hydrodynamic description of interacting fermions 
with a changing interaction parameter g x leads to a gen- 
eralized inhomogeneous Luttinger liquid actio n 17 ' 18 
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where is a canonical bosonic field. Here the effective 
velocity u x generally also changes with the interaction 
strength. We will derive the corresponding correlation 
function for an abrupt jump below. Additionally, how- 
ever, backscattering must be considered. In particular, it 
is known that the Hamiltonian density contains a leading 
relevant oscillating operator given by£ 
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where are left and right moving fermion fields and 
kp is the Fcrmi-wavenumbcr. Normally this operator 
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can be neglected under the integral, but this is no longer 
the case for inhomogeneous systems. In fact, it is ex- 
actly this operator as a local perturbation which causes 
the renormalization of defects in wires£ and spin chains £ 
The oscillating part of the interaction results in the same 
bosonic operator,— so that even a change in interaction 
alone will induce backscattcring. 

In principle conductance depends on the connections 
of the wire to both of the leads. We will concentrate here 
on the backscattering in the case where the length of the 
wire is larger than the coherence length oc u/T so that 
it suffices to consider each junction separately*^ A single 
junction can be characterized by changes taking place in 
a small region around x = 0. The field <fr x is assumed to 
be slowly varying on the scale of the Fermi wavelength 
so that it is possible to use an expansion of <f) x in Eq. ^ 
and a partial integration to derive an effectively local 
perturbation 
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while the oscillating operator cancels everywhere in the 
uniform regions. So far we have kept only the leading rel- 
evant operator, but there are higher order terms which 
will be discussed later. In general A is complex, except 
for particle hole symmetric situations. Even relatively 
smooth junctions become effectively sharper and sharper 
under renormalization, so that the relevant backscatter- 
ing is non-zero unless the amplitude A is adjusted to zero, 
which requires the fine tuning of two parameters. As we 
will see later it is indeed possible to identify such con- 
ducting fixed points in a lattice model by an appropriate 
choice of parameters on both sides of the junction. The 
existence of a conducting fixed point is also of relevance 
for the discussion about possible charge fractionalization 
in Luttingcr liquids ^ 29 ' 30 If the Luttinger liquid supports 
only charges e(l± g x )/2 related to the chiral eigenstates 
of the Hamiltonian on each side of the wire, then it seems 
to be impossible that backscattering can be tuned to zero. 
Our results thus support the analysis in Ref. [3(| that at 
such a junction an arbitrary charge can be injected into 
a Luttinger liquid. 

In order to perform a renormalization group analysis 
we have to consider the full partition function 



D(f>e 
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with the quadratic action in Eq. ([!}. For an abrupt 
change in the interaction parameter at x = from g x <o = 
9e to g x >o = g r and similarly for the velocity u x the gen- 
eral bosonic Green's function G(x,x';t) = (4> x ,t4>x' ,o)o is 
found by matching the left and right parts. In particular, 
it is possible to solve 



G(x,x';t) =T^2e iUmT G m (x,x') 
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with 
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G m (x,x') = 6(x-x') (6) 



by allowing a discontinuity in the derivative of G m at 
x = x'iH We determine the general Green's function to 
be 

G(x, x'\t) = ——In sinh 



ttT I J-L 

u x 



Ilr> 



(7) 



C[x,x']g x ^ 



sinh 


ttT (M + M - ir 




sinh 


ttT 









where we have defined a new effective interaction param- 
eter q = 2( — + — Here C\x,x'] is 1 when x and x' 
are in the same region, and when they are not. The 
renormalization of a local perturbation in Eq. ([3J can be 
determined with the help of the Green's function by inte- 
grating out the Fourier components above a cutoff A^^ 
which gives 
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We therefore expect that the effective backscattering 
renormalizcs as a power law in the temperature R oc 
T$~ l . 

As a concrete lattice model we can consider spinless 
fermions at half-filling 



H = 
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where n x = i[> x '4> x - The corresponding interaction pa- 
rameter g and the renormalized velocity u in Eq. (JT|) are 
functions of U and t, which are known from the Bethe 
ansatz.— For small jumps and interactions we have esti- 
mated the size of A perturbatively, which turns out to be 
proportional to the corresponding renormalized velocity 
field u x 
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with 2kp = 7r, so to lowest order it does not matter if the 
velocity change occurs due to inhomogeneous interactions 
or hopping amplitudes, which may even compensate each 
other. In particular, Eq. (fT0|) suggests that a conduct- 
ing fixed point can be achieved for a sharp discontinu- 
ity by a velocity matching U£ = u r , but the derivation 
above is valid only for small interactions. The numer- 
ical simulations show that this condition holds even in 
a strongly interacting model. Equation (jTU)) also implies 
that backscattering can be made arbitrarily small by very 
slow " adiabatic" changes. 

In order to calculate the backscattcring amplitude R 
numerically, we use quantum Monte Carlo simulations 3 ^ 
on systems with an abrupt junction at x = 0, i.e. U x< o = 
Ue and U x >o = U r and analogously for the hopping am- 
plitude t x . For long system sizes L > 40t x /T the bound- 
ary condition at ±L/2 becomes irrelevant. The backscat- 
tcring R induces a 2kp interference pattern in the den- 
sity, the so-called Friedel oscillationS ) 13 ' 16 i 32 which can 
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FIG. 1: (Color online) The local response in the density x for 
a jump in interaction from Ug = to U r = 1.8t r at T = 0.1t r - 
Circles (blue): No discontinuity in hopping, tg = t T . Squares 
(red): The hopping on the left is adjusted to tg ~ 1.518t r in 
order to match the velocities on both sides. Solid lines (black) 
are fits to the predicted behavior in Eq. (|14p for even and odd 
sites separately. 
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FIG. 2: (Color online) Effective backscattering amplitude 
R(T) on a logarithmic scale extracted from the amplitude 
of the density oscillations in Eq. Q14[ l for a jump from (a) 
Ue = to U r = 1.8t, (T K = 3.4 x 10~ 4 i), and (b) Ug = It to 
U r = lAt with t = tg = t r , (T K = 6 x 10~ 6 t). The amplitudes 
are extracted from fitting the local response for x < ("left 
side") and x > ("right side") separately. 



be calculated directly in the simulations and also give 
additional information about the correlation functions. 
In particular, we consider the local oscillating density in 
a half-filled lattice in response to changing the chemical 
potential 
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where the density is bosonized as 



1 (— l) x 
n x = no -= d x 4> x + const. sin[v / 47T0£ C ] (12) 



analogously to the local magnetization and susceptibility 
in spin chains 

The results for the Friedel oscillations x — Xo + 
(— l) x Xalt near an abrupt junction at x = are shown 
in Fig. Q]for a change from Ug = to U r = 1.8t r at fi- 
nite temperatures T = Q.lt r for two cases: In the case 
that the interaction strength changes but the hopping is 
equal tg = t r , we have ug < u r and strong alternating 
2kp = tt oscillations are observed (circles). In the sec- 
ond case, the hopping tg on the left was also increased 
in order to exactly match the velocities 2tg = ug = u r rj 
3.036i r (squares). Clearly, the backscattering oscillation 
is strongly suppressed with a different position depen- 
dence, but it is not zero. As the numerical data will 
show it turns out that the relevant backscattering term 
is exactly zero in this case. The uniform parts of \x 
on the two sides are constant and approximately given 
by the zero temperature result xo = ^f~j which are not 
equal in either case. We have also tested other cases 
with more complicated changes in hopping and interac- 
tions over three sites in order to tune A = in Eq. (fl"0|). 
and in all cases backscattering is strongly suppressed. 

Let us first analyze the position dependence of the 
Friedel amplitude in the case of unequal velocities, which 



shows a characteristic maximum in Fig. [T] reminiscent 
of the local susceptibility near ends in spin chains! 32 ' 35 
However, we will show that the behavior is not exactly 
of the same form as for scattering from open ends as was 
conjectured beforeji 3 - In order to calculate the alternat- 
ing response in the presence of the perturbation ([3]) with 
small |A| cx \u r — ug\ we use the bosonized form of the 
density (12]) 



f 1/T 9 I \ 
Xait = A/ dr— (sin\/47rf/). T , T cosv / 47i : (/'o,o) ,(13) 
J Oil \ /o 

where the dependence on fi is given by a shift of the 



field d x (f> by 



which turns the sine-depcndcncc 



into a cosine-dependence with an additional factor of 
x. Using the Green's function in Eq. ((7]) we arrive 
at integrals of the form J Q dr\ sinh(X — iiiT)\~ 29 = 
2 ff (sinh2X)- 9 P_ g (coth2X) where P t (z) is the Legendre 
function. Therefore we find 
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where z = coth[27rTx/u. E ]. Compared to the scattering 
from an open boundar y 13 ' 32 there are two notable differ- 
ences: first of all there is an additional factor in the form 
of the Legendre function, which changes the shape signif- 
icantly near the scatterer, but quickly approaches unity 
for x > u x /ttT. Second, there is an additional temper- 
ature dependent factor oc T^ _1 which is in agreement 
with the renormalization behavior predicted in Eq. (|8]l. 
As can be seen by the fit in Fig. [TJ the behavior (fTl)) 
describes the numerical data perfectly, where x is taken 
to be the position from an effective scattering center. We 



4 




FIG. 3: (Color online) The alternating part of the density 
oscillations at T = 0.1i r on the interacting side U = 1.8t r , 
where the hopping on the non-interacting side is adjusted so 
that ue = 0.98it r , ue = u r , and ug — 1.02u r , respectively. 



have confirmed the position dependence for many differ- 
ent temperatures and discontinuities. This unambigu- 
ously shows that the action in Eq. ([1]) and the Green's 
function (JT)) leading up to Eq. (fT4|) arc a reliable descrip- 
tion of the problem. 

An interesting detail in Fig. [T]is the maximum in the 
non-interacting region where Ue = and ge = 1 which is 
absent for an open end but now arises due to the Legendre 
function with g < 1. Physically this can be understood 
as a proximity effect, where the behavior in a range on 
the non-interacting side is influenced by the collective 
excitations on the other side. 

The temperature dependent amplitudes of the fits to 
the position-dependent part in Eq. (|T4")) directly give the 
backscattering R(T) in Fig. [5J where we have used the 
amplitude from open ends corresponding to R = 1 as 
the normalization, which is known exactly ! 35 ! 36 The data 
confirm the predicted power-law behavior at low tem- 
peratures. Note that the independent fits on both sides 
give roughly the same R(T). For the larger jump in 
the interaction U = 1.8t r there are deviations at higher 
temperatures coming from higher order operators (left 
panel). Below a characteristic temperature Tjc we ex- 
pect the power law to break down as the stable fixed 
point R(0) — 1 is approached, but this energy scale could 
not be reached in the simulations. 

In order to find a conducting scenario in the low tem- 
perature limit, it is interesting to analyze the case of 
equal velocities ui = u r with strong discontinuities in 
both hopping and interaction in more detail. As shown 
in Fig. Q] density oscillations are still observed in this 
case, but we are interested only in the contribution of 
the leading relevant operator in Eq. ([3]) , which will grow 
while the temperature is lowered and must change sign 
as a function of velocity at the conducting fixed point. 
In Fig. [3] we show the alternating part of the oscilla- 
tions for a jump from Ue = to U r = 1.8t r , where the 
hopping on the left has been adjusted to three different 



cases: ug — 0.98w r , ug — u r , and ug = 1.02w r . It is 
quite apparent that a sign change takes place exactly at 
ug = u r , which means that the relevant backscattering 
vanishes. The remaining oscillations visible in Fig. Q] for 
ug = u r are caused by higher order local operators in 
Eq. ([3j. In particular, the next- leading terms are given 
by d x e 1 ^^^ and e lv/167r< ^ with scaling dimension g + 1 
and 4<7, respectively, which are irrelevant unless Ag < 1. 
Away from half filling the marginal operator d x 4> x =o will 
also be present^ but does not affect the scattering to 
first order. The velocity matching rule for a conduct- 
ing fixed point is surprisingly simple, considering that it 
is not linked to any special symmetry in this scenario. 
Perfect conduction in quantum wires with impurities has 
been described before in cases where a renormalization 
to the periodic boundary condition fixed point occurs.— 
This, however, is not the case here. We find a novel con- 
ducting fixed point described by Eq. (J7J which does not 
correspond to a standard boundary conformally invariant 
theory. 

To conclude, we have studied the intrinsic backscat- 
tering present when connecting a quantum wire to a 
lead. All inhomogeneities and impurities in the junc- 
tion are in general relevant for repulsive Luttinger liq- 
uids, leading to a conductance which scales as G(T) ex. 
£-[1 - (T/Tk) }, with an unusual power-law exponent 
given by g = 2( + in terms of the interaction pa- 

rameters on both sides of the junction gg, g r . We have 
shown that it is possible to achieve a perfect connection, 
i.e. an absence of relevant backscattering, in inhomoge- 
neous wires by tuning other parameters such as the ve- 
locity in the lead and in the wire. The general Green's 
function was calculated in the presence of an abrupt jump 
along the wire based on an inhomogeneous Luttinger liq- 
uid action, which was in excellent agreement with numer- 
ical simulations. The results suggest that the observation 
of Friedel oscillations in the density along the wire, which 
can be attempted by scanning probe microscopy, can be 
used to analyze local scattering centers. A systematic 
study of inhomogeneities and thus an experimental test of 
our results seems feasible in semi-conductor heterostruc- 
tures. One could start from a device where a quantum 
wire adiabatically broadens into a two-dimensional elec- 
tron gas^ By applying a sharp potential barrier at the 
interface a local backscattering center could then be real- 
ized allowing it to continuously tune the backscattering 
amplitude A. 
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